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Abstract. By means of simple quantum-mechanical models we show that under 
certain conditions the main assumptions of the connected moments expansion (CMX) 
are no longer valid. In particular we consider two-level systems, the harmonic oscillator 
and the pure quartic oscillator. Although derived from such simple models, we think 
that the results of this investigation may be of utility in future applications of the 
approach to realistic problems. We show that a straightforward analysis of the CMX 
exponential parameters may provide a clear indication of the success of the approach. 

1. Introduction 

Some time ago Horn and Weinstein [1] introduced the t-expansion for the calculation of 
the ground-state energy of quantum-mechanical systems. It is a Taylor expansion about 
t = of a generating function E{t) where the coefficients are cumulants or connected 
moments of the Hamiltonian operator. The main problem posed by this approach is the 
extrapolation of the t-power series for t — )■ oo in order to obtain the ground-state energy. 
Horn and Weinstein [1] and Horn et al ^ proposed the use of Pade approximants and 
later Stubbins [3] tried some other extrapolation techniques. Without doubt, the most 
popular extrapolation strategy was proposed by Ciowsloski |:4j and it gives rise to the 
connected-moments expansion (CMX) and leads to an expression for the systematic 
calculation of the energy of the ground state. Knowles [5] studied the CMX and derived 
an elegant and compact expression for the approximants in terms of matrices built from 
the cumulants. 

It is well-known that the series of CMX approximants for the energy exhibits 
singularities and convergence problems that limit its usefulness [5HTT]. In an attempt 
to overcome those difficulties some authors proposed variants of the CMX [HHinillSHE] • 

The t-expansion, as well as the CMX and its variants, have been tested on 
several simple models with varied success [3| [9| IT2 | IT61421] . and in spite of their 
notorious limitations they have even been applied to several problems of physical 
interest pl43t[6t[9t[TT | [T3 | [T H l25 | [26]. including the calculation of the electronic energy 
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of atoms and some small molecules |^[^ [27H2U] . It is well-known that the convergence 
properties of the series of CMX approximants may be considerably poorer than those 
of the Rayleigh-Ritz variational method in the Krylov space and the Lanczos algorithm 
[5|[6|[9t [TlHT3l [23 t[24ll26] . However, they can be improved by an appropriate choice of the 
reference function [SO]. The main reason for still insisting on the development of the 
CMX appears to be its size consistency [Ij which those other approaches do not obey. 
However, in most of the applications summarized above size consistency is not an issue. 

In this paper we test the main assumptions of the CMX by means of simple 
quantum-mechanical models. In this way we expect to draw useful conclusions that 
may apply to realistic problems for which such a detailed analysis is not feasible. In 
section [2] we outline the main ideas behind this approach. In section |3]we apply the CMX 
to n-level models. In section H] we resort to the harmonic oscillator and a nontrivial 
anharmonic oscillator. Finally, in section |5] we discuss the results and draw conclusions. 

2. The Cumulant or t— expansion 

In order to facilitate present discussion, in this section we outline the main ideas 
behind the t-expansion (or cumulant expansion) and the CMX. The moment-generating 
function 

Z{t) = {^\e-'^\^) = ±t^f^, (1) 

j=o J- 

gives us the moments of the Hamiltonian operator H, fij = Iv^), in the reference or 

trial state \lp) that we assume to be normalized {(p \ip) = 1. The logarithmic derivative 
of this function 

E(t) = = (2) 

exhibits several interesting properties: 

• E{t) > Eq for all t, where Eq is the ground-state energy. 

• E'{t) < 
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• limt^oo E(t) = Eq provided that the overlap between \(p) and the ground state lipo) 
is nonzero {{ipo \(p) ^ 0). 

The function E(t) is closely related to the cumulant function K(t) defined by 
Z{t) = e^'(*) [31]. The formal Taylor series of E{t) about t = yields the t-expansion 
(also known as cumulant or cluster expansion): 



j=o J- 



(3) 



where the cumulants Ij (or connected moments) can be easily obtained from the 
recurrence relation [1] 



2 = 



li+iUj-i, j = 0, 1, 



(4) 



The main goal of the approach proposed by Horn and Weinstein [I] and Horn et al |2] is to 
find an appropriate summation method for the cumulant expansion ([3]) and extrapolate 
the resulting expression for t — )■ oo to obtain Eq. 

In order to carry out such extrapolation Cioslowsky [4j proposed an expansion in 
terms of exponentials 



E{t) = Eo + Y, Aj exp{-bjt) 



(5) 



where one obtains the adjustable parameters Eq, Aj, and bj by straightforwardly 
matching the Taylor series about t = of the left- and right- hand sides. Cioslowsky |1] 
showed that one can obtain Eq without explicitly calculating the nonlinear parameters 
bj. From the properties of the Fade approximants Knowles ^ derived an elegant, 
compact and systematic expression for the calculation of Eq: 

I ^ ^ ^ f r \ 



E, 



(m) 



'm+l 



/3 
h 



h 
h 



Im+2 
Im+3 



y Im+2 Im+3 ' ' ' -'2m+l J y J-m+1 j 

The CMX will be successful provided that limm-s>oo -E'o™'' = Eq if ('?/'o \'^) 7^ 0. 



h 
h 
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For simplicity, throughout this paper we assume that the eigenfunctions {ipj) of the 
Hamiltonian operator H 

H = E, \^,) , Eo<E,<E2<... (7) 

form a complete basis set so that 

oo 

= Y.^j iV^i) ' = (V^j W) (8) 

j=0 

provided that {ipi \ipj) = 6ij. Under such conditions we can write 
which clearly shows that lim^^oo E{t) = mm{jj{Ej, cj ^ 0}. 



We may say that the main assumption in the CMX is that equation ([5]) is valid 
provided that cq 7^ 0. In such a case the ground-state energy is given by the limit of 
the CMX approximants ([6]) for m — )■ 00. If the square matrix in Eq. ([6]) is singular, the 
corresponding approximant is omitted. 

In principle, it may happen that Z{t) = for some values of t in the complex 
t-plane. In such a case the cumulant expansion ([3]) converges for t < \ts\, where tg 
is the root of Z{t) = closest to the origin. (The location of such singular points 
will obviously depend on the trial function \(p)). Therefore, matching the Taylor series 
about if: = for the two sides of equation ([5]) may give rise to some difficulties, especially 
because the expression in the right-hand side does not take into account such singular 
points. Consequently, it is unclear that we can successfully extrapolate the approximate 
expression for E{t) thus derived to the limit t — t- 00. This aspect of the problem was 
addressed by Witte [33j and Witte and Shankar [3l] for some particular problems. In 
sections [3] and H] we will show that under certain circumstances the main assumption of 
the CMX, namely equation ([5]), does not apply to some simple problems. 

The exponential expansion ([5]) does in fact appear to take into account the singular 
points of the inverted series f32j. If we keep the first two terms, solve for t and 
differentiate with respect to E we obtain 
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This approximate result partially agrees with the exact one that we will derive for a 
two-level model in section |3l 



3. Simple models with n states 

Some time ago Knowles [5J applied the CMX to the two-level model 
H 



yV 1 J 



and concluded that the series ([6]) for the reference function \ipi) 



(A 



111 



converges 



to the exact eigenvalue for \V\ less than about 0.1, but as \V\ is increased the series 
converges to a result which is more and more in error. He also found that the CMX 
series reproduces the correct V-power series for the same reference function but the 

series can become violently oscillatory when the reference function is \ip2) = ^ 

V ^ / 

Besides, in the limit ^ — )■ alternate approximants for \(p2) tend to zero and infinity. 
Knowles [5] concluded that those results placed some doubt on the claims that the 
series is convergent for any reference function having a nonzero overlap with the exact 
wavef unction. 

Our numerical experiments for = 0.1 show that the CMX series ([6]) converges 
towards the ground-state energy Eq = —0.009901951358 when \(pi) is the reference 
function and to the excited state Ei = 1.009901951 when the trial vector is \<f2)- 
The CMX series also converges for greater values of V; for example, for = 1 we 
obtain the ground-state energy Eq = —0.6180339887 and the excited-state energy 
El = 1.618033988 with the former and latter reference vectors, respectively. The rate 
of convergence of the CMX series decreases (increases) with \V\ in the former (latter) 
case. Convergence towards the excited state is unexpected because there is no doubt 
that E{t) always tends to Eq as t ^ oo when Cq ^ 0. It may be for this reason that 
convergence of the CMX series to excited states has been ignored as far as we know. A 
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reasonable explanation for such anomalous behaviour is that the limit of the series of 
CMX approximants ([6]) is determined by the maximum overlap {{ipj \(p)\ as suggested 
by the fact that in the particular example just discussed {{ipo \ > |v^i)| and 
\ f2)\ > lii'o |v^2)| for all values of V. We will discuss this issue in more detail 
below. 

It is obvious that if the series of CMX approximants converges to an excited state 
when Co 7^ then the main assumption of the method, equation (jS]), is not valid in 
general. In what follows we explore this point in more detail and write the CMX 
approximants to E{t) as 

M 

E^''\t) = Ao,M + E A-.Me-'-"* (12) 

The two-level system is suitable for a discussion of the convergence properties of 
the t-expansion. In this case, equation ([9]) becomes 

m = (13) 

where AE = Ei — Eq and ^ = |cip/|cop. It is clear that the expansion ([3]) converges 



for all \t\ < \ts\, where 



ts = ^{ln^±m) (14) 

denotes the two singular points of E(t) closest to the origin of the complex t-plane. We 
see that the smallest radius of convergence occurs when ^ = 1. 
For large t we have the exponential expansion 

oo 

E{t) = Eo + AEu ^(-1)%^ u = ie-'^^ (15) 
and not the algebraic terms found by Witte [33] and Witte and Shankar [31] for more 



elaborate models in the thermodynamic limit. This exponential expansion converges 

for M < 1. Although the large-t expansion of E{t) exhibits the correct exponential 

behaviour, it is not valid at the matching point t = unless ^ < 1. Therefore, we expect 

that the CMX exponential expansion ( fT2l) is valid only if ^ < 1 (|co| > |ci|). 

If we solve equation f|T3|) for t and differentiate the result we obtain 
dt 1 
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that resembles the approximate expression (11 Op derived by Samaj et al [52] • By means 
of this simple model we realize why Stubbins [3] obtained reasonable results from a 
Pade analysis of the derivative dt/dE. Note that in this case t{E) exhibits logarithmic 
singularities instead of the branch cut singularity appearing in the spin-1/2 isotropic 
antiferromagnetic XY chain [33j. Besides, the derivative f lT6|) exhibits two singular points 
and the expansion in powers oi E — Eq |30] converges only for \E — Eq\ < Ei — Eq. 

In order to test the performance of the CMX we choose the even simpler problem 
given by the diagonal Hamiltonian matrix 

\ 

(17) 



H 



1 
2 



and an arbitrary reference state \ip) 



I 



\ 



(1 + 0-^/^ 

^1/2(1 + ^)-l/2 



, where < < oo. The 



CMX expansion coefficients for M = 2 are 

2^3 + 1 



^0,2 
^,2 
^,2 

&2,2 



e Re' -4e + i)ve-ioe + 1-^3 + 11^2-11^ + 1 



2(e + i)(e-e + i)(e-ioe + i) 



2 + 1) {i^ - + 12^ - lie + 1) 

ve - loe + 1 + 3 (e - 1) 

2^e + l) 

ve-ioe + i-3(e-i) 



1^ 



2(e + l) 

We appreciate that Aq^2 is closer to i?o = 1 when ^ < 1 and to i?i = 2 when ^ > 1 as 
shown by the expansions 



e<i 



^0,2 = 2-1 + 1-1 + ..., e>l (19) 

When < e < 5 — 2\/6 ~ 0.10 the nonlinear parameters bj 2 are real and positive. 
The CMX converges towards the ground-state energy and E^'^\t) exhibits the correct 
behaviour. When 5 — 2^/6 < ^ < 1 those nonlinear parameters are complex conjugate 
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each other 612 = ^22 with Re(6i.2) > 0. In this case the CMX still converges towards 
Eq and E^'^\t) is an almost reasonable approximation to E{t): it exhibits a minimum at 
some if: > and tends to a limit close to Eq ast 00. When 1 < ^ < 5 + 2-\/6 ~ 9.90 the 
nonlinear parameters are still complex conjugate each other 61^2 = ^2.2 but Re(6i_2) < 0. 
In this case the CMX converges towards the excited-state energy and E^'^^t) exhibits 
a minimum at some t < and tends to a limit close to i?i as t — t- —00. Finally, when 
^ > 5 + 2^/6 the nonlinear parameters bj^2 are real and negative. The CMX converges 
towards Ei and E^'^\t) tends to a limit close to Ei as t ^ —00. This discussion is 
illustrated in Fig. [1] that shows the real and imaginary parts of hj^2 as functions of 9, 
where ^ = < 9 < 7i/2, that is a more convenient variable for plotting the 

exponential parameters. 

As a further illustration of the discussion above in what follows we show the 
approximate function E^'^\t) in the four regions already indicated. For example, for 

e = o.i 

, 1002 50 -i3t 20 -i4t 

E^^'it) = + e— e— 

^ ' 1001 143 77 

is a monotonous decreasing function that leads to a limit close to Eq = 1 when t — )■ 00. 

When i = 0.5 



2 fVT5t\ 2^/T5 f VTEf 

- cos — - — -— sm 



-t 
e 2 



9 V 6 ; 45 V 6 
leads to 10/9 ~ 1 as t —> 00 but it exhibits a minimum at t ^ 2.8. When ^ = 1 

3 y2sin (^) 

E^'\t) = - (20) 

oscillates and does not tend to any limit. The CMX approximants Aq^m oscillate yielding 
either 3/2 or ±00 when M is even or odd, respectively. When ^ = 2 the approximate 
generating function 



2 (V^t\ 2v^ /yiSt^ 

- cos H sm 

9 I 6 / 45 I 6 



62 
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is unsuitable for t > and approaches the excited state as it tends to 17/9 ~ 2 as 

t — > — oo. Finally, when ,^ = 11 

(2).,^ ^ 2663 ( USVS _ _55_\ *(^+|) _ ( USVs _55_\ t[l-^) 
^ ' 1332 V 2664 1332 / V 2664 1332 / 

approaches the excited state as t — )■ oo as it tends to 2663/1332 2. Fig. [2] illustrates 

the behaviours of those functions. 

Table [T] clearly shows that when = 4 Aq^m already converges towards i?i = 2 as M 

increases as discussed above. We conclude that when > 1 the approximate functions 

E^^'^\t) obtained by matching the Taylor series about t = cannot exhibit the correct 

exponential behaviour for t — )• oo. It is not difficult to understand what happens when 

^ > 1 in the case of the two-level model. First of all note that \imt^-oo E{t) = Ei. 

Second, the exponential expansion of the cumulant-generating function (fT3i) 

oo 

E{t) = Ei- AEu'^ 5](-l)%-^' (21) 

j=0 

converges for all t < when > 1. Third, the match between the exponential expansion 
(IT^ and the t— expansion at t = is valid for both, positive and negative, values of 
t. Therefore, when |ci| > |co| the CMX simply chooses the convergent exponential 
expansion for t < and the result is the energy of the excited state. 

The CMX fails completely when ^ = 1 because both exponential expansions {t > 
and t < 0) are divergent at t = 0. It is worth noting that the original approach of Horn 
an Weinstein pj in terms of Fade approximants [M/M]{t) for E{t) yields accurate results 
in this most unfavourable case. These approximants exhibit either a saddle point or a 
minimum when M is either odd or even, respectively. From such stationary points we 
obtain Eq = 0.9993674107, 1.000359205, 0.9999637967, 1.000025380 with M = 5, 6, 7, 8, 
respectively. In addition to it, these Fade approximants yield increasingly accurate 
estimates of the poles of E(t) at tg = ±(2j + l)Tci, j = 0,1, . . .. 

The behaviour of the exponential parameters discussed above is not restricted to 
the case M = 2. Fig. |3] shows that two of the exponential parameters bj^s behave in a 
similar way. The third one, which we arbitrarily chose to be 63^3, is real because the 
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three parameters are solutions to a cubic equation. In fact, the exponential parameters 
bj^M are in general the M roots of the pseudo secular determinant 

I = (22) 

Fig. m shows Ao^m{G) for some values of M. We appreciate that the series of CMX 
approximants converges towards Eq = 1 when 6 < n/A and towards Ei = 2 when 
6 > Ti/A and diverges when 6 is close to 7r/4 = 1). It is a further illustration and 
confirmation of the arguments above. 

These examples clearly show that the condition cq 7^ is insufficient to guarantee 
the validity of equation ([5]) that is the main assumption of the CMX. In fact, they 
suggest that if |co| 7^ is not the largest overlap, then the CMX series may converge 
towards an excited state in which case E^^^\t) will not be an acceptable approximation 
to E{t) for t > 0, except for sufficiently small values of t. Besides, for some reference 
states (like those with ^ = 1 for the two-level model) the CMX does not converge at all. 

It is not difficult to show that the results for the two-level model also apply to 
arbitrary problems. For example, if we choose the reference state 

\ip) = c^\ip,) + c,\iPj) (23) 

then we draw conclusions similar to those above for ^ = |cjp/|cop because we already 
have a two-level problem. 

We have also carried out a numerical and analytical study of the CMX for 
Hamiltonian matrices of greater dimension. For example, if we choose a diagonal matrix 
with elements Hij = j6ij, i,j = 0,1,..., 9 and the reference vector with coefficients 
Cj = (1 + 5jo)/\/T3, j = 0,1,..., 9 then the generating function E(t) is a ratio of 
polynomial functions of u = e~*. The w-power series (large-t expansion) converges 
for all \u\ < 1.074531738 that is quite close to unity. For this reason we expect some 
difficulties when matching the t- and w-series. In fact, in this case the CMX series does 
not appear to converge towards the ground-state energy (or the convergence rate is too 
small for practical purposes). 



Further analysis of the connected moments expansion 
4. Oscillators 



12 



The conclusions of the preceding section are not restricted to the particular case in which 
only a finite number of states contribute to E(t). In what follows we consider quantum- 
mechanical systems with an infinite number of bound states. In order to keep present 
discussion as simple as possible in what follows we consider the harmonic oscillator 

(24) 

If we choose 

|(/?) = A^exp(-2a;V5) (25) 

where is a normalization factor, then the CMX converges towards the ground state 
because |co| > \cj\ for all j. In this case the exponential parameters for the case M — 3 
are 61,3 ^ 4, 62,3 ~ 8 and 63,3 12.6 that are quite close to the actual excitations 
AEj = E2j — Eq as expected. 

If, on the other hand, we choose 

\^) = n(^x'^ - exp (-2a;V5) (26) 

then the CMX converges towards the second excited state because \c2\ > \cj\. Note 
that in this case {ipo \ip) 7^ and limt^oo E(t) = 1 ; consequently E^^\t) will not be an 
acceptable approximation to E{t) beyond a sufficiently small neighbourhood of t — 0. 
This conclusion is consistent with the exponential parameters 61,3 —3.87, 62,3 ~ 4 and 
63,3 ~ 9.29. 

If we choose 

\^) = ^m + m (27) 

then the CMX series Aq^m alternate between ±00 and 3 for M odd and even, respectively. 
This result is not surprising in the light of our previous discussion of two-level models. 
At first sight one may think that it suggests that the CMX fails when two overlaps 
are equal \ci\ — \cj\ and greater than any other |cfe|; however, this is not the case. For 
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example, the reference function 

/ r040v^ 64 

^ V 4217 ^ 4217 g-.^ / 2 ^3^2 + l) _ ^2 + l) (28) 
71-4 ^ ^ ^ ^ 

satisfies cq = C2 = \/ ^113859^ ~ 113859 CMX converges towards the ground state. 

In this case we obtain 61^3 = 62.3 ~ 8-53 + 1.76i and 63,3 ~ 4.06. 
As a nontrivial example we choose the anharmonic oscillator 

+ (29) 



rf2 



dx'^ 

and the reference states 



=iVexp(-3xV2 

\^)e = iV (x^ - 1) exp (-3x72) (30) 



4. 

The CMX series converges towards the ground-state energy in the former case (6i_3 ~ 
6.34, 62,3 ~ 17.19, 63^3 ~ 34.42) and oscillates around the second excited-state energy in 
the latter one (61^3 ~ —2.61, 62,3 ~ 9.04, 63 3 ^ 26.1) as shown in Table [2j One expects 
that |cog| > \cjg\ for all j > and that \c2e\ > \cje\ for all j 7^ 2. What we know for sure 
is that E^^^\t) will not be an acceptable approximation to E{t) when \lp) = \lp)^ except 
in a neighbourhood of t = as indicated by the exponential parameters bj^^. 

5. Conclusions 

The main assumption of the CMX is that the exponential expansion ([5]) is a suitable 
approximation to the cumulant-generating function (jj]) when |co| 7^ 0. Up to now, it 
was believed that the only limitations of the method were singularities in the CMX 
approximants (E]) or an occasional divergence for some reference state. By means of 
simple quantum-mechanical models we have shown that the performance of the CMX 
is mainly determined by the overlaps \cj\. Present investigation strongly suggests that 
the series of the CMX approximants ([6]) converges towards the energy of the state with 
greatest overlap with the reference function. Thus, we have the following situations: 

• When I Co I > \cj\ for all j > the series Aq^m converges towards the ground-state 
energy and we can say that the CMX is valid. 
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• The series Aq m converges towards the energy of an excited state (or oscillates about 
it). In such a case we may say that the CMX is still useful but its main assumption 
embodied in equation (|5]) is no longer valid. This situation may happen when the 
reference state exhibits the largest overlap with an excited state. 

• The series Ao.m does not converge at all and the CMX is of no practical utility. 

In no case does the CMX converge to a meaningless limit as suggested by Knowles [5] 
based on results for a two-level model. The main conclusions of this paper agree with 
those drawn earlier that the CMX and its variants should be applied with extreme 
care in order to obtain useful results |5l-[TT|l2^. In fact, in most cases the Rayleigh- 
Ritz variational method in the Krylov space is by far preferable to the CMX and its 
variants [TT1I21] . We have also shown that the original method of Horn and Weinstein [1] 
yields accurate results in a case in which the CMX badly fails. The choice of an 
appropriate trial state may improve the results dramatically as illustrated by a recent 
study on the Rabi Hamiltonian [30] . 

Finally, one of the main contributions of this paper is that the analysis of the 

exponential parameters bj M may provide an indication on the rate of convergence of 

the CMX approximants. If all such parameters are real and positive one expects a 

reasonable convergence towards the ground state (provided that cq 7^ 0). This analysis 

is greatly facilitated by the fact that those exponential parameters are the roots of the 

relatively simple pseudo secular equation fl22|) . 
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Table 1. Convergence of the CMX towards the ground-state energy (second column) 
and excited-state energy (third column) of the two-level model ( [T7|) for ^ = 1/4 and 
^ = 4, respectively. 



IVl 


A 




2 


1.015384615 


1.984615384 


4 


1.000975609 


1.99902439 


6 


1.000061031 


1.999938968 


8 


1.000003814 


1.999996185 


10 


1.000000238 


1.999999761 


12 


1.000000000 


2.000000000 
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Table 2. Convergence of the CMX towards the ground-state (g) and second-excited 
state (e) energies of the anharmonic oscillator ([29)) 



M 


^o,m(^) 


^o,M(e) 


5 


1.060692159 


7.439371257 


10 


1.060363186 


7.456069907 


15 


1.060362073 


7.450017954 


20 


1.060362093 


7.451366303 


25 


1.060362090 


7.455118704 


30 


1.060362090 


7.454183973 


35 




7.451642486 


40 




7.454364274 


50 




7.454214745 


60 




7.453864737 


70 




7.455066766 


80 




7.455185890 


90 




7.453941990 


100 




7.453833053 


RPM pniEE] 


1.060362090 


7.455697938 
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Figure 4. Approximants Ao,m with M = 2 (line), M = 3 (squares), M = 4 (filled 
squares), M = 5 (circles) and M = 6 (filled circles) as a function of 6 



